Antiferromagnetism in a bosonic mixture of rubidium ( 87 Rb) and potassium ( 41 K) 
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We simulate the experimental possibility of observing the antiferromagnetic (AF) order in the 
bosonic mixtures of rubidium ( S7 Rb) and potassium ( 41 K) in a two-dimensional optical lattice in 
the presence of harmonic confinement. By tuning the interspecies interactions and the lattice heights 
we have found the ground states, within the mean-field approximation, that interpolate from the 
phase separation to the AF order. For a moderate lattice height the coexistence of the Mott and 
AF phase is possible for the Rb atoms while the K atoms remain in the AF-superfluid phase. This 
observation may provide an experimental feasibility to hitherto unobserved AF order for 87 Rb - 41 K 
mixture. 
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The field of ultracold atoms has set a new era for its 
outstanding achievements for simulating many quantum 
phenomena related with correlated interacting particles 
in a strong periodic potential [l[ . Such systems have been 
successfully implemented to realize nonmagnetic phases 
such as Mott, superfluid and vortex states. However, 
realizing quantum magnetism (QM) Q within current 
experimental framework is still a challenging task due to 
an extremely low value of the exchange coupling 0, Q , 
though there have been conscientious theoretical efforts 
and proposals in this direction [5l-fl(i| . The possible mile- 
stone would be the experimental realization of an anti- 
ferromagentic (AF) order in a two-component system. 

Two-component boson system with tunable inter- 
species interactions [ll[ has been a fruitful laboratory 
for studying man y in triguing phenomena rangin g fr om 
phase separation [12| and topological structures Jl3j to 
the exotic states of matter such as supersolid |9[ (su- 
perflow with broken translation symmetry) and counter- 
superfluid @. Current trends and publications for im- 
plementing such system in the regime where QM can be 
perceived are mostly theoretical, and are limited to the 
subspace of integer filling and to the hardcore bosons 
only 041] ■ Besides, the experimental realizations of ul- 
tracold atoms are always accompanied by the additional 
trapping potential that may modify the possible phase 
diagram. 

In this paper we explore the ground states of a sys- 
tem consisting of mixture of interacting bosons in a two- 
dimensional (2D) optical lattice in an external trapping 
potential. In particular, we focus on the recent experi- 
mental realizations of a degenerate mixture of rubidium 
( 87 Rb) and potassium ( 41 K) atoms [3 where one can 
exploit the inherent asymmetry of the tunneling ampli- 
tudes and interaction energies. In general, this asymme- 
try is a necessary criterion for observing the AF order 
in the two-component system 0, |i| . For a moderate lat- 
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FIG. 1: (Color online) Density plot of 87 Rb (bot- 
tom row) and 41 K (top row) for C/Rb-K = 0.28 with 
SRb = 14.0 for different number of atoms [JVjifc , JVjc] = 
[705, 110], [381, 441], [304, 548], and [288, 610] for the columns 
(a), (b), (c), and (d). The axes represent the lattice of size 
41 x 41. 



tice potential the heavier atoms (Rb) can be in the Mott 
phase while the lighter atoms (K) will remain in the su- 
perfluid phase if the interspecies interactions are weak. 
For stronger interspecies interactions, in place of a sim- 
ple phase-separation, the lattice discreteness may trigger 
magnetic order such as AF phase. 

The two-component antiferromagnetism referring here 
is the state in which the site population alternates be- 
tween two species analogous to the 2D checkerboard 
solid. By tuning the interspecies (Rb-K) interactions and 
the number of atoms of each species we have found the 
ground state where the density of each species alternate 
between adjacent sites with Rb atoms in the Mott phase 
while K atoms are in the superfluid phase. The domain 
of this checkerboard pattern extends over the size of the 
cloud when the number of atoms of each species are com- 
parable (Fig.[T](b)). For imbalanced mixtures the density 
distributions reflect the symmetry of the trapping poten- 
tial with the ring like structure of Rb atoms around the 
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densely packed K atoms as a signature of the phase sep- 
aration. 

Our formulation starts with the model Hamilto- 
nian 1151. 



^2 -^Lni l<T n i!tT > , (1) 



where the indices a ^ a' refer to Rb or K, and are 
the nearest neighbor sites. Here fit and fij i<7 are the 
creation and annihilation operators at the lattice site i 
whereas n i cr = fit a a i a is the number operator for the a 
species. The parameters J a and U a:rT ( a ,) are respectively 
the tunneling amplitudes and the onsite atom-atom in- 
teractions between a and c(cr'). The magnitude of the in- 
terspecies interactions U a ,<7' can be varied experimentally 
by using Feshbach resonance. The chemical potentials /Jo- 
are adjusted to fix the number of atoms of each species. 
The effect of the external harmonic confinement is rep- 
resented by the energy offset term in the Hamiltonian. 
In principle, Rb and K atoms may experience different 
trapping potentials due to their difference in masses. 

For the optical lattice deep enough to validate the 
model |T]) the parameters can be approximately writ- 
ten in terms of the lattice heights and atomic scatter- 
ing lengths as [r| J a ~ \As a exp(— 2.Q7^/s^)E r ^ and 
U a ^ ~ 5.97(a cr!(T /A)(s o .) - 88 i? I ., cr , where E r ^ is the recoil 
energy and s a is the lattice height of the a species. In the 
limit J a <C U a ,a an d in the absence of external trapping 
potential = 0), the Hamiltonian (TT]) can, within the 
second order perturbation theory, be mapped onto effec- 
tive spin- 1/2 Heiscnbcrg Hamiltonian In the simplest 
case when the total filling of one atom_per site the effec- 
tive Hamiltonian can be written as [11, 



(2) 



where the exchange couplings scale as t z , t± ~ J 2 jU a . 
This spin model is strictly valid when the system is in 
deep Mott regime so that the Hilbert space is spanned 



by states |n< 



T >) with n i)<T ( CT /) 



0,1. For the tunnel- 



ing amplitude close to the superfluid-Mott transition the 
model of isospin representation breaks down. However, 
this gives some hints for the parameter values where one 
should expect the magnetic order. Since the exchange 
couplings scale as t z ~ J 2 , a large tunneling would be 
beneficial for providing a large exchange coupling. 

In this study we use the decoupling mean-field (DMF) 
approximation [l7j which has been a successful model 
in obtaining the phase diagram of a single component 
system. In this approach, one can approximate the 



off-diagonal terms in the Hamiltonian as at 



(aj)a iiCT + at Ja^ a ) - (&?)(&,>) , where (a ii<r ) = fa 



is the so-called superfiuid order parameter. Within this 
approximation the mean-field Hamiltonian can be writ- 
ten as the sum of the on-site Hamiltonians Hmf = 
y^ - Hi , which can be diagonalized self-consistently us- 
ing the truncated number basis. In our simulations we 
choose the initial order parameter fa^ a as real random 
variables in the interval {0, 1}, and iterate the solu- 
tion self-consistently until it converges to a desired ac- 
curacy (10~ 8 ). In each realizations of the simulation we 
readjusted the chemical potentials for a given set of pa- 
rameters, and the number of atoms are fixed accordingly. 
We have picked the converged number of atoms in the 
range Ar^ + Nk ~ 800 to 900 when the solution relaxes 
to the ground state. In all calculations the size of the 
lattice is fixed to L 2 = 41 x 41. 

As mentioned earlier, our study concentrates on the 
mixtures of 87 Rb and 41 K in a two-dimensional optical 
lattice. For the parameters in the simulations we take 
the masses of Rb and K atoms as TORb = 87.0 u and 
tok = 41.0 u with u = 1.67 x 10~ 27 kg, and the scattering 
lengths as ajtb = 98.0 ao and ok = 63.0 ao with ao = 
0.529 x lO" 10 m. We fix the lattice height for the Rb 
atoms SRb = 14£V,Rb assuming that the lattice height 
for the K atoms can be adjusted independently. 

Figure [T] represents an example plot of the density 
of K (top row) and Rb (bottom row) atoms for a set 
of four pairs of atom numbers. We have considered 
the lattice height SRb = 14.0 for the Rb and the cor- 
responding Hubbard parameters are Ufa, = 0.30 and 
Ju,b = 0.008, whereas for the K atoms we have taken 
Uk = 0.092, and J K = 0.049 respectively. For con- 
venience, wc have taken the same trapping potentials 
for both the Rb and K ato ms, the f requency of the Rb 
atoms being uj = 2tt x 60-\/sRb/60Hz. The frequency 
for the K atoms can be found by setting the condition 
WRbW^b = ?tikw k . AH the energies are expressed in 
terms of the recoil energy of the Rb atoms. In our ob- 
servations the phase separation is likely scenario if the 
number of atoms is large and if /Tjib-K > v / ^ r Rbf^K- 
For a non-lattice system the critical value of the Rb-K 
scattering length for the onset of the phase separation 
is 2^/aRbTORbOKTOK/ (»TiRb + to k) [3- It should also be 
noted that the critical value for the superfluid-Mott insu- 
lator transition in a homogeneous system with unit filling 
is (U a /zJ a ) c = 5.8 fHf]. If we neglect the effect of the 
trapping potentials the Rb atoms remain in the Mott in- 
sulating (if the integer filling prevails) while the K atoms 
remain in the superfiuid phase. 

The density distribution in Fig. Q] form a ring of Rb 
atoms around the dense K atoms at the center of the trap. 
It is expected since Rb atoms are interacting stronger 
than K atoms and are pushed outward from the center 
of the trap. For small number of K atoms a small patch 
of AF order is formed near the center of the cloud sur- 
rounded by an extended Mott plateau of Rb atoms. As 
the number of K atoms increases the AF order extends 
over the whole lattice, and the correlator (a\a a i) shows 
minimum when A^b ~ Ak . For large number of K atoms 
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FIG. 2: (Color online) Atom distribution for the Rb (square) 
and K (triangle) atoms along the line x = 0. The data is 
extracted from the Fig. \T\ 




the overlap region shrinks in the form of ring near the 
periphery while the central region is void of Rb atoms. 
In both extreme cases we have signs of phase separation 
along with the existence of the magnetic order. 

In Fig. [5] we plot the cut of the density profile along 
the line x = for the same parameters as in Fig. [TJ The 
Mott plateau of Rb and superfluid profile of the K atoms 
exist along with the AF order in the overlapped region. 

In order to detect the hidden order in the lattice exper- 
imentally a number of proposals have been put forward. 
Two prominent techniques are the noise correlation 

and the lattice modulation spectroscopy [H[ . Here we 
emphasize the noise correlation by analyzing the images 
of the density distribution during the ballastic expansion 
of the cloud in the far field regime [2l|, [22[ . 

We define the normalized correlation function in the 
momentum space 



C(k) 



J dq(n(q)h(q + k)) 
fdq(h{q))(h(q + k)) 



(3) 



If the state of the system is known, the expectation 
values that appeared in the integrals can be found easily. 
In particular, if the state is defined by the product of the 
number state \(po) = IIj|nj) , which is a valid assumption 
for Rb atoms, the correlation function turns out to be 



C Rb (fc) = 



,Rb 



1 



N 2 

JV Rb 



N 2 

JV Rb 



nl 

fc.Rb 

M 2 

JV Rb 



(4) 



Here we have neglected the term containing the delta 
function. 

In Fig. [3] we show the density plot of the correlation 
function Cm,(k) in the k x — k y plane corresponding to 
the parameters of Fig. [1] (b). In addition to the Mott 
peaks for integer values of (k x ,k y ) the correlation func- 
tion develops additional peaks. It should be noted that 
for the checker-board solid the size of the unit cell in 



FIG. 3: Noise correlation function Citb(k) for Rb for the pa- 
rameters of Fig. [1] (b). Mott peaks of the original lattice are 
spanned by the vector k — l(nx,my) whereas the additional 
peaks are at the position k = i/v2(— nx, my) where n and m 
are integers and I — 2n/\d\, d being the lattice constant. 



the real space gets doubled and rotated by it/ 2 with re- 
spect to the original lattice. It can be easily verified that 
the length of the unit vector in the reciprocal lattice is 
l/-\/2 times that of the original lattice. Therefore, the 
secondary peaks in the correlation function do not ap- 
pear at the position half integral multiple of k x or k y , 
but at the center of the unit cell spanned by the orig- 
inal reciprocal lattice vector. For this particular set of 
parameters in Fig. [3] the ratio of the amplitudes of the 
secondary (k = 1/2) to the primary (k = 0) peaks in the 
correlation function is 0.45. 

Another possibility to observe the AF order in the two 
component system is to create an excitation by means 
of lattice shake or by modulating the intensity of the 
lattice height periodically. The AF order can be revealed 
with an additional peak of C/Rb-K in the lowest excitation 
spectrum [23[. This is because the lattice modulation 
transfer a Rb atom to the nearest occupied lattice site 
and the gain in energy will be U^b-K if the nearest site is 
occupied by the K atom. This approach can be suitable 
for a strong lattice confinement with both components 
forming Mott insulating states. But for the parameters 
we are dealing in this paper the number of K atoms in 
a given site is not well defined due to the atom number 
fluctuations and the lattice modulation may give a broad 
peak in the excitation spectrum. 

In order to refine our analysis we have also varied the 
trapping potentials for the Rb and K separately. In this 
case we have not observed a clear signature of the AF 
order but the phase separation only. This may be due 
to the fact that the difference in local chemical poten- 
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FIG. 4: (Color online) Density plot of 87 Rb (bottom row) and 
41 K (top row) for U Rh - K = 0.45 with s Rb = 20.0. The num- 
ber of atoms are [Nw,,Nk] = [188,640], [354,458], [580,200], 
and [665,100] for columns (a), (b), (c), and (d). The axes 
represent the lattice of size L 2 = 41 x 41. 

tials between Rb and K atoms greatly exceed the tiny 
exchange interaction needed to stabilize the AF order. 
In the experimental situation the equal trapping poten- 
tial may not be a problem as the frequencies can be tuned 
by setting the condition mRbW^ = ttikw^. 

To be more concrete with our speculation that the 
large hopping amplitude could be a regime conducive for 
the AF order we consider a deeper optical lattice with 
SRb = 20 for the Rb. In Fig. [4] we show the density plot 
for the typical atom numbers of each species when the in- 
terspecies interactions is C/pfb-K = 0.45. The simulation 
parameters are U^b = 0.41, Uk = 0.13, Jr^ = 0.0025, 
and Jk = 0.026. The trapping frequency for the Rb 
atoms is ui = 2ir x 60^/sRb/60Hz while the frequency of 
the K atoms is found by the equal potential conditions. 
The Mott plateau of the Rb atoms and the supcrfluid pro- 
file of the K atoms still persist but the AF order washes 
away. There are few patches of the AF order in the region 
where two species overlap. 

From the theoretical point of view when both species 
are in deep Mott regime there could be an AF order but 
the corresponding exchange couplings are extremely low. 
In the deep Mott regime atom tunneling is exponentially 



suppressed, and a tiny perturbation such as thermal and 
quantum fluctuations are sufficient to destroy the fragile 
AF order. Our observation is also consistent with the 
previous phase diagram Q of two-component bosons on 
a square lattice at half-integer filling of each species. For 
the parameters of Fig. Q] we have 2z JRb/C^Rb-K = 0.22 
and 2zJK/C^Rb-K = 1-4, which lies in the 2CB-phase of 
Fig. 1 in Ref. 

In [l4[ it is reported that the minor fraction of the K 
atoms reduces the visibility of the interference pattern of 
the Rb atoms, we speculate that this impurity-induced 
loss of coherence [3 H3] may be accounted for the pres- 
ence of the AF order, and the detail experimental findings 
in the regime of comparable atom number of each species 
would reveal more information. 

Unlike in a single component system where the chemi- 
cal potential fixes the total number of atoms, the two- 
component system allows number fluctuations for the 
same set of chemical potentials when the interspecies in- 
teractions are switched on. The asymmetry observed in 
the density distribution in Fig. [4] may be due to the de- 
generacy present in the ground state since a set of chemi- 
cal potentials may give the different set of atom numbers. 
It is interesting to compare the quantum results shown 
in Fig. U with the experimental demonstration of simul- 
taneous existence of superfluidity and magnetism in the 
spinor Bose-Einstein condensate J2j|. 

In summary, we have simulated a typical scenario of a 
degenerate bosonic mixtures of the Rb and K atoms in 
a 2D optical lattice in the presence of external trapping 
potential. By tuning the interspecies interactions we find 
that the AF order is possible if the Rb atoms are in the 
Mott phase while K atoms are in the superfluid phase. 
This ordered phase can be reflected in the density distri- 
bution as well as in the noise-correlation signal. In our 
view there has not been any study (at least to our knowl- 
edge) on the AF order in the two-component system 
when one species remain in the supcrfluid phase while 
the other in the Mott phase, this observation may pro- 
vide the novel regime for studying quantum magnetism 
in ultracold system. 
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